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Abstract 

We present a class of models for the initial mass function (IMF) for stars forming 
within molecular clouds. This class of models uses the idea that stars determine their 
own masses through the action of powerful stellar outflows. This concept allows us 
to calculate a semi-empirical mass formula (SEMF), which provides the transformation 
between initial conditions in molecular clouds and the final masses of forming stars. For a 
particular SEMF, a given distribution of initial conditions predicts a corresponding IMF. 
In this paper, we consider several different descriptions for the distribution of initial 
conditions in star forming molecular clouds. We first consider the limiting case in which 
only one physical variable - the effective sound speed - determines the initial conditions. 
In this limit, we use observed scaling laws to determine the distribution of sound speed 
and the SEMF to convert this distribution into an IMF. We next consider the opposite 
limit in which many different independent physical variables play a role in determining 
stellar masses. In this limit, the central limit theorem shows that the IMF approaches 
a log-normal form. Realistic star forming regions contain an intermediate number of 
relevant variables; we thus consider intermediate cases between the two limits. Our results 
show that this picture of star formation and the IMF naturally produces stellar mass 
distributions that are roughly consistent with observations. This paper thus provides a 
calculational framework to construct theoretical models of the IMF. 

Subject headings: stars: formation - ISM: clouds - galaxies: formation 



1. INTRODUCTION 



The initial mass function (IMF) is perhaps the most important resuh of the star 
formation process. A detailed knowledge of the initial mass function is required to un- 
derstand galaxy formation, the chemical evolution of galaxies, and the structure of the 
interstellar medium. Unfortunately, however, the current theory of star formation says 
very little about the IMF (see, e.g., the reviews of Shu, Adams, & Lizano 1987, hereafter 
SAL; Zinnecker, McCaughrean, & Wilking 1993). In particular, we remain unable to 
calculate the initial mass function from first principles. 

Given the extreme importance of the IMF and the many successes of the current 
theory of star formation, we feel that it is now time to begin building models of the IMF. 
The purpose of this present paper is to present a class of IMF models which use the idea 
that stars, in part, determine their own masses through the action of powerful stellar 
winds and outflows (see, e.g., SAL; Lada & Shu 1990). Within the context of the current 
theory of star formation described below (§1.2), we can conceptually divide the process 
which determines the IMF into two subprocesses: 

[1] The spectrum of initial conditions produced by molecular clouds (the star forming 
environment). 

[2] The transformation between a given set of initial conditions and the properties of 
the final (formed) star. This transformation is accomplished through the action of 
stellar winds and outflows. 

Notice that molecular clouds are not observed to be collapsing as a whole; on average, the 
lifetime of a molecular cloud is (at least) an order of magnitude longer than the free-fall 
time (e.g., Zuckerman & Palmer 1974). Thus, these clouds exhibit quasi-static behavior 
and it makes sense to conceptually divide the process of determining the distribution of 
stellar masses into the two steps given above (see Zinnecker 1989, 1990). 

A large body of previous work on the IMF exists in the literature (see, e.g., the 
reviews of Zinnecker, McCaughrean, & Wilking 1993; Elmegreen 1985). Many of these 
studies use the idea that fragmentation of clouds leads directly to the masses of the 
forming stars (e.g., Hoyle 1953; Larson 1973; Bodenheimer 1978; Elmegreen & Mathieu 
1983). More recent work (Larson 1992, 1995) has extended these ideas to include the 
observed fractal and hierarchical structure of molecular clouds (e.g., Scalo 1985; Dickman, 
Horvath, & Margulis 1990; Scalo 1990; Lada, Bally, & Stark 1991; Houlahan & Scalo 
1992). Zinnecker (1984, 1985, 1989, 1990) has discussed the two subprocesses given 
above and has explored several different fragmentation schemes to produce the IMF. The 
concept that stars help determine their own masses has just now begun to be incorporated 
into models of the IMF. Silk (1995) has discussed the IMF for stars which have masses 
limited by feedback due to both ionization and protostellar outflows. Nakano, Hasegawa, 
& Norman (1995) have introduced a model in which stellar masses are sometimes limited 
by the mass scales of the formative medium and are sometimes limited by feedback. 
Finally, a more primitive version of this current theory has been presented previously 
(Adams 1995). 

For the point of view of the IMF adopted in this paper, traditional arguments based 
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on the Jeans mass are not applicable. A characteristic feature of molecular clouds is 
that they are highly non-uniform; dumpiness and structure exist on all resolvable spatial 
scales. In particular, no characteristic density exists for these clouds and hence no (single) 
Jeans mass exists. We stress that, at least in the context of present day star formation in 
molecular clouds, the Jeans mass has virtually nothing to do with the masses of forming 
stars. 

1.1 The IMF Observed 

We begin this discussion by emphasizing that stars can only exist in a finite range 
of masses. Stellar objects with masses less than about 0.08 Mq cannot produce central 
temperatures hot enough for the fusion of hydrogen to take place; objects with masses 
less than this hydrogen burning limit are brown dwarfs (see, e.g.. Burrows, Hubbard, & 
Lunine 1989; Burrows et al. 1993; Laughlin & Bodenheimer 1993). On the other end 
of the possible mass range, stars with masses greater than about 100 Mq cannot exist 
because they are unstable (e.g., Phillips 1994). Thus, stars are confined to the mass range 

0.08<m<100, (1.1) 

where we have defined m = M^/{1Mq). Notice that this mass range is rather narrow in 
the sense that it is much smaller than the conceivable range of masses. Stars form within 
galaxies which have masses of about 10^ ^Mq and stars are made up of hydrogen atoms 
which have masses of about 10~^^ g ~ 10~^^Mo. Thus, galaxies could build objects 
anywhere in the mass range from lO^^^M© to I^^^^Mq., a factor of 10^^ in mass scale. 
And yet, as we have discussed above, stars live in the above mass range which allows 
stellar masses to vary by only a factor of ~ 10^. 

The initial mass function in our galaxy has been estimated empirically. The first 
such determination (Salpeter 1955) showed that the number of stars with masses in the 
range m to m + dm is given by the power-law relation 

f (m) dm ^ dm , (1-2) 

where the index b = 2.35 for stars in the mass range 0.4 < m < 10. However, more 
recent work (e.g.. Miller & Scalo 1979; Scalo 1986; Rana 1991; Tinney 1995) suggests 
that the mass distribution deviates from a pure power-law. The distribution becomes 
flatter (and may even turn over) at the lowest stellar masses {b approaches unity for the 
lowest masses 0.1 < m < 0.5) and becomes steeper at the highest stellar masses {b ~ 
3.3 for m > 10). The observed IMF can be approximated with an analytic fit using a 
log-normal form (Miller & Scalo 1979), i.e., 

logio /(logio m) = ao- ai log^g m - a2(logio m)^ , (1-3) 

where qq = 1.53, ai = 0.96, and 02 = 0.47. The true IMF has more structure than 
a simple log-normal form (Scalo 1986; Rana 1991), although equation [1.3] provides a 
good analytic reference distribution. Figure 1 shows three successive approximations to 
the observed IMF: the Salpeter power-law [1.2], the Miller/Scalo log-normal form [1.3], 
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and the more recent distribution taken from Table 2 of Rana (1991). We note that 
the construction of the IMF from observational quantities (e.g., the observed luminosity 
function) requires considerable processing. However, the basic features of the IMF seem 
to be very robust. As a general rule, the IMF does not change very much from one star 
forming region to another. For the sake of definiteness, in this paper, we use the analytic 
fit given by equation [1.3] as a benchmark with which to compare our theoretical models. 

1.2 The Current Theory of Star Formation 

In the last decade, a generally successful working paradigm of star formation has 
emerged (see, e.g., SAL for a review). Since the IMF models of this paper use this 
paradigm as a starting point, in this section we quickly review its basic features. One 
result of this present work is thus a consistency check - we show that this star formation 
paradigm can produce an IMF similar to that observed. 

In our galaxy today, star formation takes place in molecular clouds. These clouds 
thus provide the initial conditions for the star forming process. Molecular clouds have 
very complicated substructure. In addition, these clouds exhibit molecular linewidths 
Av which contain a substantial non-thermal component (e.g., Myers & Fuller 1992); this 
linewidth broadening is generally interpreted as a "turbulent" contribution to the velocity 
field. 

Molecular clouds are supported against their self-gravity by both "turbulent" mo- 
tions and by magnetic fields. The fields gradually diffuse outward (relative to the mass) 
and small centrally condensed structures known as molecular cloud cores are formed. 
These cores represent the initial conditions for protostellar collapse. In the simplest pic- 
ture, these cores can be (roughly) characterized by two physical variables: the effective 
sound speed aes and the rotation rate O. The effective sound speed generally contains 
contributions from both magnetic fields and "turbulence", as well as the usual thermal 
contribution. The total effective sound speed can thus be written 

2 2 2 2 / \ 

^eff = ^therm + ^mag + ^turb • (1-4) 

The molecular cloud cores eventually undergo dynamic collapse, which proceeds 
from inside-out; in other words, the central parts of the core fall in first and successive 
outer layers follow as pressure support is lost from below (Shu 1977). Since the infalling 
material contains angular momentum (the initial state is rotating), not all of the infalling 
material reaches the stellar surface. The material with higher specific angular momentum 
collects in a circumstellar disk. The collapse flow is characterized by a well defined mass 
infall rate M, the rate at which the central object (the forming star/disk system) gains 
mass from the infalling core. Notice that no mass scale appears in the problem, only a 
mass infall rate. In particular, the total amount of mass available to a forming star is 
generally much larger than the final mass of the star. 

One important characteristic of the rotating infalling flow described above is that 
the ram pressure of the infall is weakest at the rotational poles of the object. The central 
star/disk system gains mass until it is able to generate a powerful stellar wind which 
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breaks through the infall at the rotational poles and thereby leads to a bipolar outflow 
configuration. Although the mechanism which generates these winds remains under study 
(see the review of Konigl & Ruden 1993), the characteristics of outflow sources have been 
well studied observationally (see the review of Lada 1985). One of the basic working 
hypotheses of star formation theory is that these outflows help separate nearly formed 
stars from the infalling envelope and thereby determine, in part, the flnal masses of the 
stars (SAL; Lada & Shu 1990). In this paper, we use this idea as the basis for calculating 
a transformation between the initial conditions in a molecular clouds core and the flnal 
mass of the star produced by its collapse (see §2; Shu, Lizano, & Adams 1987, hereafter 
SLA; Adams 1995). 

1.3 Organization of the Paper 

This paper is organized as follows. In §2, we derive a semi-empirical mass formula 
which provides a transformation between the initial conditions in a star forming region 
and the flnal masses of the stars formed. In subsequent sections, we use this transforma- 
tion in conjunction with the observed properties of molecular clouds to derive an initial 
mass function. In §3, we flrst consider the limit in which only one physical variable (the 
effective sound speed) determines stellar masses. In §4, we consider the opposite limit in 
which a large number n of physical variables contribute to the determination of stellar 
masses; in the limit n ^ cxd, the central limit theorem implies that the IMF approaches 
a log-normal form. In §5, we explore more complicated models of the IMF; these mod- 
els are intermediate between the limiting cases studied in the two previous sections. In 
this section we also consider the effects of binary companions on the IMF. Finally, we 
conclude in §6 with a summary and a discussion of our results. 

2. A SEMI-EMPIRICAL MASS FORMULA 

In this section, we calculate the transformation between initial conditions and the 
flnal masses of the stars produced. In other words, we derive a semi-empirical mass 
formula (SEMF) for the masses of forming stars. In the current picture of star formation, 
the flnal masses of stars are produced in part through the action of powerful outflows. 
Thus, we must know how this stellar outflow stops the inflow and separates the star/disk 
system from its molecular environment. Although this process has not been well studied, 
we obtain a working estimate by balancing the ram pressure of the outflow against that 
of the infall. For this calculation, we adopt the arguments flrst presented by SLA (see 
also Adams 1995). 

The key concept in this argument is that the flnal mass of a star is determined by 
the condition that the stellar outflow is strong enough to reverse the direct infall onto 
the star. We write this condition in the form 

M„ = (2.1) 

where is the mass loss rate of the wind and M* mass infall rate onto the star itself; 
this infall rate is generally only a fraction of the total mass infall rate because much of 
the infalling material falls directly onto the disk. Notice that we should really compare 
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the ram pressure (~ M^v^) of the wind with that of the infall (~ M*v*). However, both 
velocities are determined by the depth of the steUar potential well and are thus compa- 
rable in magnitude; we thus divide out the velocities and incorporate any uncertainties 
into the parameter 6. 

Since we do not yet fully understand how high velocity outflows are produced, we 
must proceed in a semi-empirical manner (although considerable progress in this area has 
recently been made - see Shu et al. 1988, 1994). The kinetic energy i?out of the outflow 
will generally be some fraction a of the binding energy of the star, i.e., 

Eout = . (2.2) 

-ft* 

The natural time scale associated with stellar processes is the Kelvin-Helmholtz time 
scale. We thus take the duration of the outflow (which is produced by a stellar process) 
to be a fraction /? of the Kelvin-Helmholtz time, i.e.. 

Tout = P-B-r ■ (2.3) 

Combining the above two equations, we thus reproduce the observational correlation that 
the mechanical outflow luminosity Lout = -E'out/T'out is roughly a constant fraction of the 
photon luminosity of the central source, i.e., 

a 

Lout ~ 'p ^*'> (2-4) 

where observations show that a//? ~ 10~^ and that this correlation holds over several 
decades of (see Bally & Lada 1983; Lada 1985; Levreault 1988; Edwards, Ray, & 
Mundt 1993). Finally, if the winds roughly conserve energy while driving bipolar outflows, 
we obtain the result 

M^-—^e-L„ (2.5) 
-ft* p 

where e is an additional efficiency parameter. 

The strength of the infall can be measured by the rate M* at which matter falls 
directly onto the star. The total infall rate onto the central star/disk system is given 
by the collapse solution for an isothermal cloud core (Shu 1977). This infall rate M for 
purely spherical infall takes the form 

M = moa^/G, (2.6) 

where mo = 0.975 is a dimensionless constant. For cloud cores which are not isothermal, 
the mass infall rate can still be written in the form [2.6] with the sound speed a taken to 
be the total effective sound speed (see equation [1.4]) and a different numerical constant 
(Adams et al. 1995). 

When rotation is present, not all of the material falls all the way in to the stellar 
surface (see Cassen & Moosman 1981; Terebey, Shu, & Cassen 1984). The material with 
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higher specific angular momentum collects in a circumstellar disk whose radius is roughly- 
given by the centrifugal radius 

^o^^^. (2.7) 

When the stellar radius is small compared to the centrifugal radius, i?* <C Rc, the direct 
infall rate M* onto the star itself is a small fraction of the total mass infall rate and is 
given by 

M. = = ^7^f'' . (2.8) 

The first equality is taken from equation [24] of Adams & Shu (1986); the second equality 
follows from the expressions for M and Re- 
in general, the stellar mass is only a fraction of the total mass M that has 
collapsed to the central star/disk system at a given time; we write this condition in the 
form 

=7M. (2.9) 

Disk stability considerations (Adams, Ruden, & Shu 1989; Shu et al. 1990) greatly limit 
the allowed range of the fraction 7. We generally expect 7 ~ 2/3; smaller values of 7 
imply larger relative disk masses and hence systems that are gravitationally unstable. 

The combination of the above results implies that the final properties of the newly 
formed star are given by the following SEMF: 

L,M! = S,n,,Hl^ = A^,. (2.10) 

where we have defined a new dimcnsionless parameter A in the second equality. Under 
most circumstances, we expect that the parameters a, /3, 7, 5, and e can be estimated 
to a reasonable degree of accuracy. For example, disk stability arguments suggest that 7 
~ 2/3 and empirical estimates imply that /?/q! ~ 10^. The parameter e is close to unity, 
whereas the parameter 5 has a value of a few. We thus expect that the parameter A will lie 
in the range 10^ < A < 10^. In real star forming environments, the parameters a, /3, . . . 
do not have exactly the same values for all cases. Instead, the values of these parameters 
have a distribution which is determined by the underlying physics of the problem. 

Equation [2.10] provides us with a transformation between initial conditions (the 
sound speed a and the rotation rate Q) and the final properties of the star (the luminosity 
and the mass M*). If we use "typical" values for present day clouds (e.g., a = 0.35 
km and ~ 3 x 10~^^ rad s~^ ~ 1 km s"-*^ pc~^) and the observed protostellar 
luminosities (L* ~ 2OL0), we obtain stellar mass estimates ~ 1 Mq, which is the 
typical mass of stars forming in regions with these properties. In spite of its highly 
idealized nature, the SEMF [2.10] thus provides reasonable estimates for the mass M* 
as a function of initial conditions {a,Q). It is useful to write this transformation in 
dimcnsionless form ^ 

Lm^ = 20A3 45l)^^ (2.11) 

where we have defined L = L^/{1Lq), m = M^/{1Mq), 035 = a/(0.35 km s"-*^), Qi = 
km s"^ pc"^), and finally A3 = A/IO^. 
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Notice that much of the uncertainty in this calculation has been encapsulated in 
the parameter A, which should really be considered as a complicated function of all the 
environmental parameters. Notice also that we have set the final mass of the star by the 
criterion that its outflow is sufficiently powerful to overwhelm the infall; in actuality, the 
star will continue to gain some mass, both from residual infall and from disk accretion, 
after this evolutionary state has been reached. This mass correction can be absorbed 
into the factor A in equation [2.10]. We have also characterized the initial conditions by 
only two physical variables (a and O) whereas much more complicated initial states are 
possible. Finally, we have ignored radiation pressure in this argument; for sufficiently 
massive stars (M* > 7 Mq) radiation pressure will help the outflow reverse the infall 
(e.g., see Wolfire & CassineUi 1986, 1987; Nakano 1989; Jijina & Adams 1995). 

In order to evaluate the semi-empirical mass formula derived above, we must deter- 
mine the relationship between mass and luminosity for young stellar objects. In general, 
the luminosity has many contributions (Stabler, Shu, & Taam 1980; Adams & Shu 1986; 
Adams 1990; Palla & Stabler 1990, 1992). For our present purposes, however, we can 
simplify the picture considerably. The most important source of luminosity for low mass 
objects is ultimately from infall; in other words, infalling material falls through the grav- 
itational potential well of the star (and disk) and converts energy into photons. This 
luminosity can be written 



where 035 and m are the dimensionless sound speed and mass as deflned above. For the 
stellar radius, we have used the scaling relation i?* = (3 x 10^^ cm ) 035, indicated by the 
stellar structure calculations of Stabler, Shu, and Taam (1980). The efficiency parameter 
r] is the fraction of the total available energy that is converted into photons. For spherical 
infall, all of the material reaches the stellar surface and r] ^ 1. For infall which includes 
rotation, some of the energy is stored in the form of rotational and gravitational potential 
energy in the circumstellar disk. We generally expect ~ 1/2. 

In addition, the star can generate its own internal luminosity through deuterium 
burning, gravitational contraction, and eventually hydrogen burning. This additional 
luminosity contribution is important for stars with masses larger than ~ few Mq. For 
low mass stars on the main sequence, the luminosity is a very sensitive function of stellar 
mass, L* ~ Mf. For higher mass stars, the mass/luminosity relation flattens to the 
form ~ (and eventually flattens further to the form ~ M* for very high mass 
stars). For stars still gaining mass from infall, the internal luminosity contribution is 
somewhat different, but has been calculated for much of the relevant range of parameter 
space (Stabler 1983, 1988; PaUa & Stabler 1990, 1992; see also Fletcher & Stabler 1994). 
For our present purposes, we use the following simple approximation for the internal 
luminosity 



LA = r] 



GMM 



70Lq iial^m, 



(2.12) 




(2.13) 



which is roughly valid for the mass range IMq < M* < IOM0. 
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Putting both contributions together, we obtain the luminosity as a function of mass 
in dimensionless form: 

L = lQr]al^m . (2.14) 

Thus, at low masses, the luminosity is dominated by the contribution from infall and L 
is a linear function of m. At higher masses L ~ with the crossover point at m ~ 3.3 
(for representative values of 77 = 1/2 and 035 = 1). At sufficiently high masses, the form 
[2.14] is no longer valid and L ~ lOOm^ for the range 10 < m < 100. 

Finally, putting all of the above results together, we present the SEMF in dimen- 
sionless form: 

m^Qm[Kz/riY'^ al^il~'^^^ low m , (2.15a) 
m = 1.65 Ag^^ 035^^ f];^ intermediate m, (2.15b) 

m = Qm afJ^VL'^'^ high m. (2.15c) 

This SEMF provides a transformation between initial conditions and the final mass 
of the star. One way to view this result is shown in Figure 2. Here, we assume that the 
effective sound speed 035 and the rotation rate VLi are the two most important parameters 
which determine the initial conditions. We thus set all of the remaining parameters to 
constant values such that A3 = 1 and we take r] = 1/2. Figure 2 shows the resulting 
contours of constant mass in the plane of initial conditions, i.e., the (035, Oi) plane. The 
region in the far upper left corner of the diagram corresponds to brown dwarfs, i.e., 
objects with masses less than the hydrogen burning limit. The region in the lower right 
corner corresponds to stars that are too massive to be stable; the lower right part of the 
diagram also corresponds to initial conditions for which radiation pressure helps limit the 
stellar mass (see Jijina & Adams 1995). In actual star forming regions, the parameters 
(in addition to 035 and Oi) which enter into the SEMF will have a distribution of values 
(roughly centered on the values assumed here). As a result, the set of initial conditions 
which lead to a star of a given mass will be a band in the (035, plane instead of a line. 
Finally, for comparison, we note that previous authors have considered the mean density 
n and the mass infall rate M as the two most important variables which determine stellar 
masses (see Figures 1 - 3 of Nakano et al. 1995). 

3. EMPIRICAL MODEL: THE INITIAL MASS FUNCTION 
FOR CLUMPY MOLECULAR CLOUDS 

In this section, we consider the limiting case in which the effective sound speed is 
the only physical variable which determines stellar masses. Here, we use two observed 
scaling laws to determine the distribution of the effective sound speed and hence the 
distribution of initial conditions for star formation. This result, in conjunction with the 
SEMF derived in the previous section, produces a nearly power-law IMF in reasonable 
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agreement with observations. We also calculate the efficiency of star formation from this 
model. 



3. 1 Observed Distributions of Initial Conditions 

Observational work spanning many different star forming regions and many different 
size scales suggests that the effective sound speed in molecular cloud cores (or clumps) 
obeys a simple scaling law (see, e.g., Larson 1981; Scalo 1987; Myers & Fuller 1992). 
For sufficiently large size scales (r ~1 pc) and low density n < 10^ cm~^, the observed 
linewidths in cores have a substantial nonthermal component which scales with density 
according to the law 

Avcxp"^/^. (3.1) 

Although this result was obtained from observational data, relations of this type can be 
calculated theoretically from the supposition that magnetohydro dynamic waves (e.g., 
Alfven and magnetoacoustic waves) are the source of the non-thermal motions (see 
Fatuzzo & Adams 1993; McKee & Zweibel 1995). Notice that this scaling law is valid over 
a finite range of densities; at sufficiently large densities the observed linewidths become 
equal to the thermal linewidths (in other words, the total linewidth does not vanish 
as p ^ oo as implied by equation [3.1]). 

If this velocity Av is interpreted as a transport speed, then a "turbulent" or "non- 
thermal" component to the pressure can be derived (Lizano & Shu 1989; Myers & Fuller 
1992) and has the form P = PqIyv^p/ po). This "turbulent" equation of state also implies 
a scaling relation between the line- width Af and the mass M^i of the clump. Using 
hydrostatic equilibrium arguments, we obtain the relation 

Mci ~ {^vY , (3.2a) 

where ~ 4 for the law given by equation [3.1]. Since we consider the linewidth Au to 
define the effective sound speed of the region, we can write this scaling relation in the 
form 

Mci = Mais > (3-2b) 

where the mass scale M ~ 7 Mq is determined by the normalization of the observed 
scaling law (see, e.g., Larson 1981; Myers & Fuller 1992). 

Given the above result, we must now determine the distribution of clump masses 
Mci. Many groups have studied the observed clump mass spectrum of molecular clouds 
and have found nearly power-law forms, i.e., 

where the index of the power-law typically has the value p ^ 3/2 (e.g., see Scalo 1985; 
Lada, Bally, & Stark 1991; Blitz 1993; Tatematsu et al. 1993). This scaling relation must 
have a cutoff at both high mass (to keep the total mass of the cloud finite) and at low 
mass (to keep the total number of clumps finite). 
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Although the distribution [3.3] was obtained from observations, such distributions 
can, in principle, be calculated theoretically. For example, simple models which envisage 
clouds to be composed of a collection of interacting clumps (Norman et al. 1995) can 
be used to derive clump mass distributions of the general form [3.3]. We emphasize that 
much more theoretical work on this subject should be done. 

3.2 A Simple Model for the IMF 

We can now piece together all of the above arguments to construct an initial mass 
function. Here, we interpret the line-width Av as the effective transport speed a which 
determines the initial condition for star formation. We use the semi-empirical mass 
formula of §2 to provide the transformation between the initial conditions and the final 
stellar properties. To start, we consider the simplest case in which only the sound speed 
varies and the remaining parameters of the SEMF are kept constant. The relationship 
[3.2] between clump mass and linewidth, in conjunction with the clump mass spectrum 
of equation [3.3], determines the distribution of the effective sound speed. Combining 
this distribution with the SEMF, we obtain an initial mass function of the form 

diV _ dN dMci 
^ ~ dM^ ~ dMci dM^ 

where we have assumed that only a single star forms in a given clump. The power-law 
index b of the distribution is given by 

b = q{p-l)/ii+1^2/ii+l, (3.5) 

where /U is the scaling exponent which determines how the stellar mass varies with effective 
sound speed. As shown by equation [2.15], this index lies in the range 11/6 < /x < 
11/3. Thus, this simple argument produces a power-law IMF with an index in the range 
6 = 1.6 — 2.1. This result compares reasonably well with the observed power-law index 
of the IMF which has b ^ 2.35 (see Salpeter 1955). 

In Figure 3, we show the IMF calculated from this model. Here we use equations 
[3.2] and [3.3] to determine the distribution of the effective sound speed. We then use 
the SEMF in the form of equation [2.11] and the mass/luminosity relationship [2.14]. 
The result is shown as the dashed curve in Figure 3; as indicated by equation [3.4], 
this distribution has a power-law index b 1.6 at low masses and b 2.1 at higher 
masses. Also shown for comparison is the analytic fit to the observed IMF (from Miller 
& Scalo 1979). Notice that the agreement between the theory and the observations is 
quite reasonable, but is not exact. We interpret this finding to mean that this picture of 
the IMF is basically correct, but it still incomplete. 

The basic logic of this model can be summarized as follows. Molecular clouds produce 
a distribution of initial conditions for star formation. In the simplest picture considered 
here, the clouds produce a distribution of clump masses. Because larger clumps have 
larger effective sound speeds due to turbulence and other small-scale physical processes, 
this distribution of clump masses implies a corresponding distribution of effective sound 
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(3.4) 



speeds, which represent the initial conditions for star formation. We then use the idea 
that outflows help determine the final masses of forming stars to find a transformation 
between the initial conditions and the final stellar masses. Using both this transformation 
and the set of initial conditions, we thereby obtain the IMF. 

3. 3 Efficiency of Star Formation 

We can directly calculate the efficiency of star formation from this model of the IMF. 
Here, we define the star formation efficiency £ to be the ratio of the mass in stars to the 
total cloud mass, i.e.. 



= / 77^^*^^-'/ / :7^^cic^Mei, (3.6) 



where Mi and M2 are the lower and upper cutoffs of the clump mass distribution [3.3]. 
Using the SEMF of the previous section (for simplicity, we use only the low mass version 
[2.15a]) in conjunction with the scaling law of equation [3.2], we can evaluate this integral 
to obtain the efficiency 

£ = 0.66 [A^/v]'/' fh-^'^ m^^/^-' , (3.7) 

1 + 2,/q-p 

where we have assumed that Mi ^ M2. We have also defined rh = M/{1Mq) and 
7712 = M2/(1M0). If we use representative values, A3 = 1, ry = 1/2, p — 3/2, q — 4, in = 
7, and m2 = 1000 (e.g., see Williams, de Geus, & Blitz 1994), we obtain a star formation 
efficiency £ ~ 0.07. 

The star formation efficiency calculated here must be compared with the observed 
value for giant molecular clouds taken as a whole. This efficiency typically has a value 
of a few percent (e.g., Duerr, Imhoff, & Lada 1982; Lada, Strom, & Myers 1993). We 
conclude that this semi-empirical model produces a star formation efficiency in reasonable 
agreement with observations. 

For completeness, we calculate the efficiency of star formation for an individual 
clump. This efficiency £^ in the low mass regime is then given by 

^* = = ^[^3/r/]^/-^a^5^ - 0.12a3t , (3.8) 

where the final approximate equality was obtained using the representative values A3 = 1, 
T] = 1/2, q = 4, and m = 7. A similar formula can be derived for the regimes of 
intermediate and high mass stars. 

In this model, we have assumed that only a single star forms within a given molecular 
clump. In many cases (see Zinnecker et al. 1993; Lada et al. 1993; Hillenbrand et al. 
1993), a cluster of stars forms within a single clump and the star formation efficiency of 
the individual clump can be much higher, f^ciust ~ 0.2 — 0.5. This present model does 
not take cluster formation into account, although this issue is important and should be 
addressed in future studies. 
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4. RANDOM MODEL: THE INITIAL MASS FUNCTION 
AS A RESULT OF THE CENTRAL LIMIT THEOREM 



In this section, we consider the Umit in which a large number of physical variables 
is required to determine stellar masses. We thus adopt a statistical approach to the 
calculation of the IMF. We start with the semi-empirical mass formula of §2 and consider 
it to be a product of random variables. In the limit that the number n of random variables 
is large {n oo) and the variables are completely independent, the IMF approaches a 
log-normal distribution. This result is a direct consequence of the central limit theorem 
(see, e.g., Richtmyer 1978; Parzen 1960). For the more realistic case of a finite number 
of not-completely-independent variables, we must study how the resulting distribution 
differs from a log-normal distribution (see also §5). The idea of using the central limit 
theorem to obtain a log- normal distribution has been discussed previously in models 
where the stellar masses are determined by fragmentation (Larson 1973; Elmegreen & 
Mathieu 1983; Zinnecker 1984, 1985). In this paper, we adopt a different approach using 
the SEMF of §2 as the starting point for our calculation. We also note that the number n 
of physical variables is actually finite and hence departures of the IMF from a log-normal 
form are important. 

The semi-empirical mass formula can be written in the general form of a product of 
variables 

n 

M. =n«., (4.1) 

where the aj represent the various quantities on the right hand side of equation [2.10], 
i.e., a, P, 7, S, e, a, and Q (taken to the appropriate powers). In this present discussion, 
we regard these quantities aj as a collection of n random variables. Thus, by taking the 
logarithm of this equation, we find that the logarithm of the mass is a sum of random 
variables, 

n 

In = In aj + constant , (4.2) 

where the constant term includes all the quantities in the SEMF that are truly constant 
(e.g., the gravitational constant G). 

The stellar mass is thus determined by a composite random variable that is given by 
the sum of random variables. The central limit theorem shows that the distribution for 
the composite variable always approaches a normal (gaussian) distribution as the number 
n of variables approaches infinity. In order to use this result, we must redefine the basic 
variables Inaj so that the new variables have zero mean, i.e., 

CjIMMj^O, (4.3) 

where fj is the probability density of the jth variable. The distribution fj is, in general, 
not a normal (gaussian) distribution. These new variables are related to the old 
variables aj through the relation 

= Inaj — (Inoij) = ln[Q;j/Q:j] , (4-4) 
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where angular brackets represent averages. Keep in mind that the averages are taken 
over the logarithms of the aj and not over the variables aj themselves, i.e., 

/oo 
In aj fj (In aj ) cZ In aj . (4-5) 
-oo 

Similarly, the distributions are the distributions of Inaj and not the distributions of 
aj. Each of the rescaled variables ^j has a variance aj given by 



oo 



- / ei/.(0)^0- (4.6) 

•/ — oo 

Next we construct a composite random variable (, defined by 

n n 

c^E^^- =Ew«.-M]- (4.7) 

i=i i=i 

In terms of this new variable, the semi-empirical mass formula becomes 

= Mc e^ (4.8) 

where Mc is a characteristic mass scale. The distribution of stellar masses (the IMF) is 
thus determined by the distribution of the composite variable C. The mass scale Mc is 
determined by the mean values of the logarithms of the original variables aj, i.e., 

n n 

Mc = n exp[(lna,)] = J] ' (4-9) 

where we have defined aj = exp[(lnQ;j)]. 

As we find below, the variance (a) of the composite variable ( essentially determines 
the width of the stellar mass distribution and is thus of fundamental importance for this 
present discussion. The variance is given by 

^ ^ POO POO 

(^)' = E E / / ^0^6 a /(e,-, a) , (4.10) 

j=lk=l-^-°°'^-°° 

where fiCj^Ck) is the joint probability distribution of the variables ^j and ^k- If the 
variables are statistically independent, then the joint probability is the product of the 
individual probability distributions and the integral in equation [4.10] can be separated. 
Using the fact that each of the variables ^j has zero mean (equation [4.3]) and the 
definition [4.6], we thus obtain the total variance, 

n 
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We can now define a new random variable 



C-^, (4.12) 

which has zero mean and unit variance. The central limit theorem tell us that the 
distribution of the composite random variable ( approaches a normal distribution in the 
limit n — > oo, i.e., 

/(C)^A^e-^"/^ (4.13) 

where J\f — l/\/2n is a normalization constant. This result is independent of the initial 
distributions fj. 

The mass of the star (from equation [4.1]) is related to this new variable ( through 
the relation 

InM* = lnMc + ((T)C, (4.14) 

where the variable ( now has a known (gaussian) distribution. Combining equations 
[4.13] and [4.14], we can write the distribution / of stellar masses in the form 

ln/(lnm) - ^^{ln[m/mc]}', (4.15) 

where A is a constant and where we have defined m = M^/{1Mq) and mc = Mc/{1Mq). 
Notice that the constant A just sets the overall normalization of the distribution. The 
shape of the distribution is thus completely determined by the mass scale mc and the 
total variance (a). 

In this limit, the distribution of stellar masses has exactly the same form as the 
Miller/Scalo approximation of equation [1.3]. The three constants ai, a2, and as of 
the Miller/Scalo law are related to the physical variables A, mc, and (a) through the 
relations 



opinio 

^ = aolnlO+^- f«4.65, (4.16a) 

4a2 



{af ^2.45, (4.16b) 

2a2 



ai In 10 ^ „^ / . . ^ V 

Inmc = ~ -2.35 . (4.16c) 

2a2 

Thus, the "characteristic mass scale" of the distribution is mc ~ 0.095. 

In this limit, where the SEMF involves a large number of statistically independent 
variables, we obtain a "pure" log-normal distribution. Since the observed IMF can be 
approximately fit by a log- normal distribution (see Figure 1), this model of the IMF is in 
reasonable agreement with observations. In this limit, the only relevant parameters are 
the total width of the distribution (determined by {a)) and the center of the distribution 
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(determined by mc)- The product of the mean values of all the of relevant variables 
combine to determine the mass scale mc of the distribution (see equation [4.9]). Similarly, 
the widths of all of the original variables combine to determine the total width {a) of the 
final distribution (see equation [4.11]). 

We thus obtain an important consistency check on this model of the IMF: The 
total width (cr) and the characteristic mass scale mc can be calculated and compared 
with the values required to fit the observed IMF (see equation [4.16]). The quantities 
(cr) and mc are determined by the distributions of all of the physical variables in the 
problem. In a complete theory, we could calculate these initial distributions from a priori 
considerations. In the absence of a complete theory, however, we can use observations 
of the physical variables to estimate their distributions and hence calculate (cr) and mc- 
Such a calculation is performed in Appendix B. Using estimates of the distributions 
of the observed physical variables aefj, fi, etc., we obtain values (cr) 1.8 and mc ~ 
0.25. Although these values are somewhat higher than the values required to fit the 
Miller/Scalo IMF (see equation [4.16]), we consider the agreement to be quite good, 
given the crudeness of the calculation. 



5. INTERMEDIATE EXAMPLES AND APPLICATIONS 



The previous discussion has considered the two limiting cases in which the number 
n of physical variables that determine stellar masses is either one (§3) or infinite (§4). In 
the first case, we obtain a nearly power-law IMF; in the second case, we obtain a log- 
normal IMF. In realistic star forming regions, however, we must consider the intermediate 
cases with 1 < n < oo. In particular, we must determine the form of the composite 
distribution (the IMF) for these intermediate cases. In general, the answer depends on 
the distributions of the initial variables. In this section, we explicitly calculate theoretical 
IMFs for several different cases with various distributions of initial conditions. We also 
consider the problem of binary companions. Although the theory discussed thus far only 
applies directly to the formation of single stars, we can show (§5.3) that the inclusion of 
binary companions does not greatly change the resulting IMF. 



5.1 Uniform Distributions 



In this subsection, we consider the simple case of n fundamental variables aj, each 
with the same distribution fj{\naj). We also consider the simplest type of distribution 
in which = Inaj is uniformly distributed in an interval [—w,w], i.e., 

= otherwise . 



The variance aj of each individual variable is related to the width w of the interval by 
the expression 

a] = . (5.2) 
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For a given number n of variables, we can thus obtain the required total width (a) of the 
distribution by taking 

w' = l{a)\ (5.3) 

where we use equation [4.16b] to set the value of (a). We also use equation [4.16c] to set 
mc and hence the center of the distribution. 

The resulting distribution (IMF) is shown in Figure 4 for the case n — 10. To obtain 
this result, we have used a random number generator to produce the SEMF variables 
(distributed according to equation [5.1]) and have calculated a million (10^) realizations 
of the mass. Notice that these random variables do not have a gaussian distribution, but 
the sum of random variables comes rather close to a log-normal distribution and hence 
reproduces the Miller/Scalo IMF quite well. In other words, n = 10 is "close enough to 
oo" for the central limit theorem to apply and hence for the distribution of ( to be nearly 
gaussian, i.e., for the IMF to have nearly a log- normal form. 

Next, we would like to determine how the composite distribution changes with the 
number n of fundamental variables ^j. Although the result depends on the initial distri- 
butions of the variables, we can get some feeling for this problem by using the uniform 
distribution [5.1] and varying the number n. The result is shown in Figure 5 for the 
cases n=l, 2, 3, and 5. The mass scale mc has been set to correspond to that of the 
Miller/Scalo IMF (shown as the solid curve). Notice that the composite distribution 
converges toward the log-normal limit rather rapidly with increasing values of n. As 
expected, the largest departures are for high masses, i.e., for the tail of the distribution. 



5.2 Power-law Distributions 



In this subsection, we explicitly consider the case in which all of the physical variables 
appearing in the SEMF have power-law distributions. This case is expected to be a 
reasonable approximation for many astrophysical systems. 

We consider each fundamental variable to have a distribution of the form 

^=^"r. (5.4) 

where C is the normalization constant. If the power-law index p > 1, as we assume here, 
then we must specify the lowest value of the parameter aoj, i.e., we introduce a lower 
cutoff but not an upper cutoff. Next, we write the distribution in terms of the logarithmic 
variable ^j, which is defined by 

= (P - 1) ln(Q!j/Q!oi) - 1 • (5-5) 
The corresponding distribution fj is then given by 

/,=e-(«^+i) , (5.6) 

i.e., we obtain an exponential distribution. Notice that this new variable has zero 
mean and unit variance. 
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For this case, we can explicitly calculate the composite distribution, and hence the 
IMF, from the initial distributions. Here we employ standard methods from probability 
theory (e.g., Richtmyer 1978; Parzen 1960). The Fourier transform Xj of the distribution 
fj {xj is generally known as the characteristic function of the variable is defined by 
the usual integral 

^dee-^^V(^+^) (5.7) 

where we have evaluated the integral to obtain the second equality. 
Next, we define a composite variable C according to 



1 



n . 

which also has zero mean and unit variance. The Fourier transform xc of the composite 
distribution is then the product of the Fourier transforms of the individual distributions 
evaluated at \/\/n. Thus, the composite transform xc is given by 

-in r p^-^/ -| n 

X<(A) = [x,(A/V5)] = [j^-^] . (5.9) 

To obtain the composite distribution itself fn{C)i we simply take the inverse Fourier 
transform of equation [5.9]. We thus obtain 



/n(C) = (;^^{^+ ^C}''"'e-("+-'<) , (5.10) 

i.e., we obtain a gamma distribution for the composite variable (. It is straightforward 
to show that in the limit n — > oo, the distribution [5.10] approaches a log- normal form 
(as the central limit theorem requires). However, for the special case of power-law initial 
distributions, we have an exact form for the distribution for intermediate cases (i.e., for 
finite values of n) . 

The total variance (a) of the composite distribution depends on the power-law indices 
p and the number of variables through the relation 

W^=(^. (8.11) 



Similarly, the central mass scale mc of the distribution is given by 

n 

n v^, 
— + _^lnao, 

?" = 1 



mc = + y Inctoj. (5-12) 



In Figure 6, we show the composite distributions for various values of n. Also 
shown for comparison is the log- normal form which corresponds to the limit n ^ oo. All 
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distributions shown here have the same total variance (cr) and characteristic mass scale 
mc as given by equation [4.16]. Notice that the convergence to the log-normal form is 
much slower than for the case of uniform distributions considered in the previous section. 
In particular, for sufficiently large masses, the IMF falls off like a decaying exponential 
in the variable In m (a power-law in the variable m) instead of like a gaussian. We note 
that the observed IMF seems to have a power-law tail at high masses, although the exact 
slope is somewhat uncertain (see, e.g., Massey et al. 1995 for a good discussion of this 
issue). 



In this section, we address the issue of binary companions. Thus far, the discussion 
of this paper has focused on the formation of a single star. On the other hand, most 
stars live within binary systems (see, e.g., Abt & Levy 1976; Abt 1983; Duquennoy & 
Mayor 1991; Bodenheimer, Ruzmaikina, & Mathieu 1993). Several different mechanisms 
for producing binary companions have been proposed, including capture via star/disk 
interactions (Clarke & Pringle 1991), formation from gravitational instabilities in cir- 
cumstellar disks (Adams, Ruden, & Shu 1989), and fragmentation during protostellar 
collapse (Boss 1992; Bonnell & Bastien 1992). Unfortunately, however, a complete the- 
ory of binary formation has not yet been obtained. As a result, we must once again 
proceed in a semi-empirical manner. 

The SEMF of §2 can be interpreted as providing the final mass Mp* of the primary 
star in a binary system. If we can write this SEMF in the form of equation [4.1], then 
the mass of the secondary Ms* can also be written in the general form of a product of 
variables, i.e.. 



where as = Ms*/Mp^ is the mass ratio of the binary system. Note that 0:5 < 1 by 
definition. A complete theory of binary formation would give us a theoretical estimate 
of the distribution of the mass ratio 0:5. As mentioned above, however, we do not yet 
have a complete theory of binary formation. As a result, we must use observations to 
measure and/or constrain the distribution of the mass ratio. For a given distribution of 
asi we can determine the consequences for the IMF using the framework developed in 
this paper. 

We first note that if the distribution of masses for the primary star has a known 
form (e.g., the log- normal form of §4), then the distribution of the secondary masses has 
nearly the same form. This result follows directly from equation [5.13] for the SEMF 
and holds for any as distribution that is not overly pathological (this statement can be 
made mathematically more precise - see Richtmyer 1978; Parzen 1960). Furthermore, 
the width {(j) g of this secondary distribution is directly determined from the width {a) p 
of the primary distribution through the relation 



5.3 The Effects of Binary Companions 



n 



Ms* = as n 



(5.13) 



{a)l = {cr)% + als, 



(5.14) 
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where aaS is the variance of the mass ratio distribution. In obtaining this result, we have 
assumed that the distribution of the mass ratio is independent of the other variables 
in the problem. If the ratio is not completely independent, then the expression [5.14] 
provides an upper limit to the variance of the secondary distribution. Notice that the 
width of the secondary distribution is always wider than that of the primary distribution. 
Similarly, the characteristic mass scale mcs of the secondary distribution is given by 

mcs/mcp = e<i""«>, (5.15) 



where mcp is the mass scale of the primary distribution and where angular brackets 
denote averages over the distribution. 

As a reference point, we consider the simplest case in which the mass ratio as is 
uniformly (randomly) distributed over the interval [0,1]. For this case, (TqiS — 1 and 
(IncKg) = — 1. Although the observed distribution of mass ratios is not completely uni- 
form, these values for aas and (In 05) represent reasonable estimates (see, e.g., Duquen- 
noy & Mayor 1991 for a more detailed discussion). 

The total mass distribution, including both primary and secondary stars, is the sum 
of the two individual distributions. For the sake of definiteness, we use the large n limit 
for which both distributions obtain a nearly log-normal form. In this case, we obtain the 
total distribution in the form 

/ = AA{e-^V2 ^ j,^-ic+Co?/2B- 1 ^ (5_;L6) 

where M is the normalization constant, is the binary fraction, and C, — Inm/mcp is the 
composite variable for the primary mass distribution (see equations [4.7] and [4.8]). We 
have also defined a parameter B which represents the ratio of the widths of the secondary 
and primary mass distributions, 

B^[l + als/{<Tfpf\ (5.17) 
and a parameter (q which determines the difference in the centers of the two distributions, 

Co^-^. (8.18) 

In the limit S — > 1 and (0 — > 0, we recover the primary distribution for the IMF. For 
relatively small departures from this limit, the joint distribution / is nearly the same as 
the primary distribution. For the reference case of a uniform distribution of the mass 
ratio ctg, we obtain values B ^ 1.2 and Co ~ 0.64. 

We can quantify the difference between the joint distribution / and the original 
distribution /o for the masses of the primary stars. We define an error functional Eji, 

/oo 
\fo-f\'dC, (5.19) 
-00 
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where we have normahzed the integral such that = 0] = 1. The size of the error 

estimate Er is thus a measure of how the distribution / differs from the original distri- 
bution /o- For the joint distribution given by equation [5.16], we evaluate this functional 
to obtain 

Er = + - 2^2(1 + e-CoV2(i+B^) | . (5.20) 

Using the uniform distribution to estimate Co and B, we find Er 0.04. We thus 
conclude that the efi^ects of binary companions on the IMF are not overly large for the 
paradigm considered in this paper. 

In some sense, the error estimate obtained above is overly conservative because it 
includes the differences between the two distributions over the entire mass range < m < 
00. If we normalize the two distributions to unity for a mass of 1.0 Mq and only consider 
the expected mass range for stars, the difference is much smaller. This result is shown in 
Figure 7. The solid curve shows the primary distribution /o with width and mass scale 
consistent with the Miller/Scalo estimate. The dashed curve shows the effect of adding a 
distribution of binary companions according to equation [5.16] with Co = 0.64, B = 1.2, 
and binary fraction = 0.75. 

6. SUMMARY AND DISCUSSION 

6. 1 Summary of Results 

In this paper, we have presented models of the initial mass function using the idea 
that stars, in part, determine their own masses through the action of stellar winds and 
outflows. Our results can be summarized as follows: 

[1] We have presented a semi-empirical mass formula (SEMF) for forming stars (see 
equation [2.10]). This result determines the transformation between the initial con- 
ditions for star formation and the final masses of forming stars and uses the idea that 
stars determine their own masses through the action of stellar winds and outfiows 
(Figure 2; see also SLA). 

[2] We have presented an empirical model for the IMF in the limit that the effective 
sound speed is the most important physical variable which determines stellar masses. 
In this limit, the spectrum of initial conditions for star formation is given by the com- 
bination of the observed clump mass distribution [3.3] and the relationships between 
clump mass, density, and linewidth (equations [3.1] and [3.2]). This distribution in 
conjunction with the SEMF produces a nearly power-law distribution of masses of 
forming stars (see Figure 3). This theoretical distribution is in reasonable agreement 
with the observed IMF. 

[3] The empirical model of the IMF also allows us to estimate the overall star formation 
efficiency 8 ~ 0.07. This calculated efficiency (see equation [3.7]) is in reasonable 
agreement with observations. 
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[4] We have studied the IMF in the hmit where a large number n of physical parameters 
play a role in determining stellar masses. The central limit theorem shows that for 
any SEMF which can be written as the product of parameters (as in equation [4.1]) 
the resulting distribution of stellar masses (the IMF) approaches a log-normal dis- 
tribution as the number of parameters n — > cxd. Since the observed IMF is crudely 
given by a log-normal distribution, this theory is in reasonable agreement with ob- 
servations as long as the number n of parameters which characterize star forming 
environments is sufficiently large. 

[5] When the central limit theorem applies (item [4]), the resulting IMF is specified 
by two quantities: the total width (a) of the distribution and the characteristic 
mass scale mc- For a given SEMF, both of these quantities can be calculated from 
the original parameters in the problem. The total width (cr) is determined by the 
quadrature sum of the variances of the distributions of all of the input parameters of 
the problem (see equation [4.11]). The mass scale mc is determined by the average 
of the logarithms of the input parameters (see equation [4.9]). For the SEMF of 
§2, the values of (cr) and mc estimated from observed distributions of the input 
parameters are in basic agreement with those required to fit the observed IMF (see 
Appendix B). 

[6] We have studied the IMF resulting from the physically realistic case of intermediate 
numbers of variables (1 < n < oo). In this case, the exact form for the IMF 
depends on the distributions of the original physical parameters of the problem. 
When these input parameters have uniform (flat) distributions, the convergence of 
the IMF to a log-normal form is quite rapid (Figures 4 and 5). For the case of 
power-law distributions of the initial variables, the convergence is much slower and 
the IMF retains a power-law tail at high masses (Figure 6). 

[7] We have briefly considered the effects of binary companions on the theory of the 
IMF presented in this paper. We show that the inclusion of binaries does not greatly 
change the resulting IMF; in addition, the manner in which binaries change the IMF 
can be directly calculated (see §5.3 and Figure 7). 

[8] The combination of all of these results demonstrates the consistency of the hypothesis 
that winds and outflows help determine the masses of forming stars by limiting the 
infall. 

6.2 Discussion 

At this point, we must carefully assess what we have calculated and what we have 
not. We have studied the idea that stars, in part, determine their own masses through 
the action of strong stellar winds and outflows. The question thus becomes: Have we 
proven this conjecture or have we presented mere speculation? The results of this paper 
represent much more than the latter, but unfortunately much less than the former. We 
have demonstrated the plausibility of our hypothesis, but we have certainly fallen short of 
a full proof. In this work, we have demonstrated that the idea of stars determining their 
own masses through the action of stellar outflows is in fact compatible with the observed 
distribution of stellar masses. 
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Perhaps the most important result of this work is that it provides a framework 
to approach the calculation of the IMF. In this framework, the calculation of the IMF 
involves the two steps outlined in the introduction: (1) the selection of initial conditions 
and (2) the transformation between a given set of initial conditions and the final mass 
of the star (see also Zinnecker 1989, 1990). Although treatments of both of these steps 
have been given here, the calculation of each of these steps can be refined considerably 
(see §6.5). 

The semi-empirical model (§3) and the random model (§4), as presented here, are 
opposite limits of the same underlying problem. In the empirical limit, the distribution 
of sound speed is the distribution of initial conditions. In the opposite limit of a random 
model with a large number of independent variables, the total distribution of initial con- 
ditions is independent of the distribution of the sound speed except for its contribution 
to the overall width of the distribution. The actual physical case lies between these two 
extremes, i.e., the IMF should be determined by seweraZ variables (say, n = 3 — 10) which 
are not completely independent. However, even in the limit of a single variable - the ef- 
fective sound speed - we obtain a nearly-power-law distribution which is reasonably close 
to the observed IMF. Furthermore, relatively few independent variables are necessary to 
"round out" the distribution to be even closer to the observed distribution (see Figure 
5). 

Next, we must consider the issue of uniqueness. In the limit that a large number 
of physical variables play a role in the star formation process, the central limit theorem 
implies that the resulting composite distribution (the IMF) always approaches a log- 
normal form. As a result, many different theories can, in principle, predict very nearly 
the same IMF. The best way to discriminate between competing theories is thus to look 
for the deviations of the theoretically predicted IMF from a pure log-normal form. As a 
general rule, the IMF will deviate most from log-normal at the tails of the distribution, 
i.e., at the low- mass and high- mass ends. It is thus crucial to obtain tight observational 
constraints on the IMF at both high and low masses. Unfortunately, however, the IMF 
is notoriously difficult to determine at both the high mass end (Massey et al. 1995) and 
the low mass end (Tinney 1995). This issue represents a challenge for the future. 

6.3 Implications of the Theory 

The picture of the IMF promoted in this paper can be tested, or at least highly con- 
strained, by observations. In addition, this theory makes several preliminary predictions, 
which we discuss below. 

The first issue is that of the semi-empirical mass formula. This result provides 
a transformation between the initial conditions for star formation and the final stellar 
properties. In this picture, the final stellar mass is most sensitive to the total effective 
sound speed and the relation has the form ~ with 2 < < 3. This theoretical 
result is in good agreement with the observed correlation between the mass M^^ax of the 
largest star in a region and the observed line-width (Aw) in that region (Myers & Fuller 
1993); the observed relation can be written 
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and is valid for the mass range 0.1 < m < 30. The theoretical relation M* ~ a'^ is a 
direct result of the hypothesis that stars help determine their masses through the action 
of stellar winds and outflows. Furthermore, the SEMF has the same general form for a 
variety of cases and hence this result is fairly robust (see Appendix A). 

We have shown that as the number of variables in the SEMF increases, the form of 
the IMF approaches a log- normal distribution, provided only that many different variables 
play a role in the SEMF. The degree to which an exact log-normal distribution is obtained 
is illustrated by Figures 3-7. Thus, one prediction of this theory is that the IMF should 
have (nearly) a log- normal form for any star forming environment. Furthermore, when 
the IMF deviates from a pure log-normal distribution, it is expected to have a power-law 
tail at high masses as illustrated by Figure 6. 

The mass distribution can be characterized by two parameters: the mass scale mc 
and the total width (a) of the distribution. For a given SEMF, the values of these 
parameters (cr) and mc can be calculated from the distributions of the original variables in 
the problem (see equations [4.9] and [4.11]). Although the underlying distributions might 
not be known exactly, the width and central values of the distributions may be estimated. 
Such estimates, in conjunction with the results of this paper, may be useful in determining 
the IMF for models of galaxy formation, cooling flows, and other astrophysical systems. 

6.4 The Low-Mass End of the IMF and Brown Dwarfs 

Another important issue is the lower mass cutoff for stars. The search for brown 
dwarfs has interested astronomers for many decades, both as a limiting case of stellar 
evolution (e.g.. Burrows & Liebert 1993; Laughlin & Bodenheimer 1993) and as a source 
of dark matter in the galactic halo (e.g., Heygi & Olive 1989; Adams & Walker 1990; 
Salpeter 1992; Graff & Freese 1995). In addition, microlensing experiments are now 
providing an important probe of these low mass stellar populations (e.g., Alcock et al. 
1993; Aubourg et al. 1993). Within the paradigm of star formation invoked here, however, 
the formation of large numbers of brown dwarfs is difficult. In the following discussion, 
we examine this statement in more detail for both limits presented in §3 and §4. 

In the limit that the effective sound speed is the most important physical variable 
(§3), the SEMF implies that in order to form very low mass stars, the effective sound 
speed must be very small. However, even the lowest temperatures expected in present 
day molecular clouds T ~ 10 K lead to the formation of stars with masses greater than 
the brown dwarf limit. Roughly speaking, stars need to become reasonably large in order 
to produce winds sufficiently powerful to reverse the infall; objects with masses smaller 
than the brown dwarf limit can only form within star forming regions with very small 
mass infall rates (such a scenario has been advocated by Lenzuni, Chernoff, & Salpeter 
1992; see also Zinnecker 1995). We thus expect brown dwarfs to be rare. 

In the opposite limit in which many different physical variables conspire to produce 
a nearly log-normal distribution (§4), the characteristic mass scale and total width must 
have given values {mc ~ 0.1 and (cr) ^ 1.6) in order to be consistent with the observed 
IMF. This result implies that the number of stars with masses less than mc (and hence 
objects below the brown dwarf limit) is very highly suppressed. We stress that this claim 
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is stronger than a blind extrapolation of the observed IMF into the unknown: In the 
limit of large n, the distribution approaches a log-normal form and, other than nic and 
(a), there are no additional parameters to specify. For the IMF of equation [4.14], the 
fraction of the total mass that resides in stars with masses less than the brown dwarf 
limit (taken here to he ruBD — 0.08) is ~ 5%. Notice also that this putative brown 
dwarf population corresponds to the low mass end of the usual stellar population and not 
the halo population. If brown dwarfs make up a substantial fraction of the mass of the 
galactic halo, then they must arise from a population with an IMF different from that of 
field stars. 

6.5 Future Work 

Many directions for future research along these lines remain. In this paper, we have 
presented a basic framework which can be used to calculate theoretical models of the 
IMF. Thus far, we have considered only extremely simple models for both the SEMF and 
the distribution of initial conditions. Thus, essentially all steps of the calculation can be 
improved significantly. 

For the SEMF, more elaborate models of the type considered here can be derived 
(see also Appendix A). In addition, future work should use hydrodynamic simulations to 
study the manner in which protostellar outflows reverse the infall. Such work will help 
determine or constrain the form of the SEMF. Finally, one can also work backwards from 
the observed distributions of initial conditions and see what type of SEMF is required to 
produce the observed IMF. 

For the distributions of initial conditions, both theoretical and observational ap- 
proaches should be pursued. Almost all of the relevant physical variables appearing in 
the SEMF can be measured in actual star forming regions. The corresponding distri- 
butions can also be determined, e.g., the effective sound speed a (see §2), the rotation 
rate Q (Goodman et al. 1993), the wind efficiency parameters P/ae (Lada 1985), and 
the star/disk mass fraction 7 (Beckwith & Sargent 1993; Beckwith et al. 1990; Adams, 
Emerson, & Fuller 1990). Although some data on these distributions currently exist, 
much more is necessary to fully understand the problem. In addition to finding more 
accurate descriptions of these distributions, future observational surveys can also deter- 
mine how the distributions of the variables change from one star forming environment to 
another. 

The relative importance of the physical variables is given by the size of the contribu- 
tion (Tj to the total width of the distribution (see Table 1). The effective sound speed is 
thus the most important variable. The rotation rate Q, the wind efficiency factors (3/ae, 
and star/disk mass fraction 7 are the next most important. Future studies should pri- 
oritize their efforts accordingly, both for observational approaches (as described above), 
and theoretical studies, which we discuss next. 

The distributions of the physical variables can also be calculated theoretically. For 
example, the distribution of effective sound speed can be calculated from the combination 
of the line- width vs density relationship (equation [3.1]) and the clump mass distribution 
(equation [3.3]). These distributions, in turn, can be calculated from MHD wave consid- 
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erations (Fatuzzo & Adams 1993; McKee & Zweibel 1995) and molecular cloud models 
(Norman et al. 1995). Similarly, theoretical models of the protostellar wind mechanism 
(Shu et al. 1994) and disk stability calculations (Laughlin 1994) will eventually determine 
the distributions of the dimensionless parameters a, 7, and e appearing in the SEMF. 
Calculations of this type remain in their infancy; much more work must be done in order 
to fully understand these distributions. 

In summary, we have presented a basic calculational framework which can be used 
to build theoretical models of the IMF. This approach is based on a SEMF and the 
underlying distributions of the physical variables which enter into the star formation 
problem. The simplest cases of these models are presented here and show reasonable 
agreement with the observed IMF. In the future, all steps of this calculation can be 
improved and we hope to eventually obtain a fundamental understanding of the IMF. 
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APPENDIX A: ALTERNATE MODEL FOR THE 
SEMI-EMPIRICAL MASS FORMULA 



In this appendix, we discuss an alternate derivation of the semi-empirical mass for- 
mula of §2. We show that the final result is very similar to that derived in the text. We 
thus conclude that the general form of the SEMF is fairly robust. 

At a given time in the collapse, most of the newly falling material falls to radii near 
the centrifugal radius Rc- We therefore consider the condition that the ram pressure of 
the outflow is sufficiently strong to reverse the infall at the radius Rc- This condition 
defines the effective end of the infall phase and can be written in the form 

M^v^ = SMvc , (Al) 

where vc ~ (GM/Rc)^/'^ is the infall speed at the centrifugal radius and where 5 is a 
dimensionless parameter. 

We use this condition in place of equation [2.1] (the old defining equation for the 
end of infall) and keep the remaining assumptions concerning the scalings of outflow 
strengths, etc. After some rearrangement, we obtain the alternate scaling relation 

where we have defined the parameter A ~ 500. Once again, the stellar properties appear 
on the left hand side of the equation and the initial conditions appear on the right hand 
side. 

In dimensionless form, this transformation can be written 

L R^/^ = 25OOA3 al^n:[^ , (A3) 

where we have introduced R = R^/{1Rq) and A3 = A/10^. 

Although this transformation appears to have a somewhat different form than that 
derived in §2 in the text, it leads to exactly the same scalings for the case of (low mass) 
young stellar objects with luminosity dominated by infall energy. In this case, the SEMF 
takes the form 

m = Aal^n~^^^ , (A4) 

where the parameter A contains all of the original parameters of the problem. This form 
is exactly the same as that obtained in §2. 
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APPENDIX B: CALCULATION OF THE 
EXPECTED WIDTH OF A LOG-NORMAL IMF 



In this Appendix, wc estimate the total width (a) and the mass scale mc of the 
IMF based on the SEMF of §2 and observed distributions of the fundamental physical 
parameters. 

The most important physical variable is the effective sound speed, which has the 
largest contribution to the total width. The distribution of the effective sound speed can 
be determined from the clump mass distribution [3.3] and the relation [3.2] between the 
effective sound speed and the clump mass. The result is 

^=AAa-[(^-i)''+il, (Bl) 
da 

where A/" is a normalization constant and g pa 4 and p ^ 3/2 are the indices of the 
distributions [3.2] and [3.3]. Since the index appearing in the distribution [Bl] is ~ 3, we 
need to introduce a lower cutoff ao to keep the distribution bounded. We take ao ~ 0.20 
km/s, corresponding to the thermal sound speed at a temperature T = 10 K. Next, we 
let a; = Ina/ao and calculate the mean of the variable x: 

POO -1 

(x) = (p-l)q / xdxe-^P-^^'^'' = . — . (B2) 

Jo ip - 1)9 

Thus, the relevant reduced variable which shows how the sound speed enters into the 
IMF is given by 

where is the exponent of the sound speed appearing in the SEMF (/x = 2 - 3, depending 
on the mass range). Another straightforward integration then gives us the variance 

where we have used /x = 2.5, p — 1.5, and q — 4 to obtain the numerical estimate. Thus, 
the sound speed contributes a little over half of the total width of the distribution (recall 
that {af = 2.45). 

Similarly, we can calculate the appropriate mean value of the sound speed from this 
distribution. As described in the text, we must calculate the average of the logarithm of 
the variable, and then exponentiate the result. We thus obtain 

a = exp[(lna)] = oq exp[l/(g(p - 1))] 0.33 km s"^ . (B5) 



Of the remaining variables, the dimensionless factors a, /?, 7, etc. appearing in the 
SEMF can all be treated the same approximate manner. Suppose a given variable Zj 
varies by a factor J^j and enters into the SEMF with an exponent /ij. Then, the quantity 
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= Hjlnzj lies in the interval [— yuj In^^, InjFj]. The contribution to the total 
variance is then approximately given by 

a] = (B6) 

In the following table, we list estimates for the exponents iij and the expected vari- 
ance factors Tj for the variables in the SEMF. The final columns show the contribution 
cr| to the total width of the distribution and the mean value (as defined by equation 
[4.5]) which determines the characteristic mass scale. For the exponents /Uj, we use the 
the form [2.15a] for the SEMF. The distribution for the rotation rate O was derived from 
the results of Goodman et al. (1993). The distribution for the mechanical luminosity 
factors (P/ae) for protostellar outfiows was taken from Figure 7 of Lada (1985); here 
we consider the three wind parameters a, P, and e to be coupled (not independent) and 
hence described by a single distribution. The star /disk mass fraction 7 has been esti- 
mated from various calculations of the stability of self-gravitating star/disk systems (see, 
e.g., Adams, Ruden, & Shu 1989; Shu et al. 1990; Laughlin 1994; Woodward, Tohline, 
& Hashisu 1994). Finally, the efficiency factor rj for protostellar luminosities has been 
estimated from calculations of protostellar structure (e.g.. Stabler, Shu, & Taam 1980; 
Palla & Stabler 1990, 1992) and from radiative transfer models of protostellar spectral 
energy distributions (Adams & Shu 1986; Adams, Lada, & Shu 1987; Kenyon, Calvet, & 
Hartmann 1993). 

Using the values given in Table 1 and the SEMF, we can calculate the total width 
and characteristic mass scale. We find that the calculated width is (cr) pa 1.81, which is 
slightly higher than the observed value of (cr) = 1.57. Similarly, the calculated mass scale 
is mc ~ 0.25, which is again higher than the observed value of mc = 0.095. Thus, the 
calculated values are approximately correct, but still differ from the observed values by 
a significant amount. In any case, this calculation is meant to be illustrative rather than 
definitive. 



Table 1. Mean and Variance of Fundamental Parameters 



parameter 


fij 






a 


035 


5/2 


N.A. 


1.56 


0.94 


Qi 


2/3 


3 


0.54 


1.5 


P/ae 


1/3 


10 


0.60 


100 


7 


1 


2 


0.48 


0.5 


5 


1/3 


2 


0.05 


1.0 


V 


1/3 


2 


0.05 


0.5 
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FIGURE CAPTIONS 

Figure 1. Observed estimates of tbe initial mass function. Dasbed curve sbows 
tbe power-law IMF of Salpeter (1955). Solid curve sbows tbe log-normal analytic fit 
to tbe IMF from Miller & Scalo (1979). Tbe dasbed curve witb symbols sbows a more 
recent empirical estimate of tbe IMF taken from Rana (1991). All tbree distributions are 
normalized to unity at m = 1 (M* = IMq). 
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Figure 2. The masses of forming stars as a function of initial conditions for the SEMF 
of §2. This figure shows contours of constant mass in the plane of initial conditions. The 
effective sound speed 035 = a/(0.35 km s~^) constitutes the horizontal axis and the 
rotation rate Oi = 0/(1 km s~^ pc~^) constitutes the vertical axis. The remaining 
parameters of the problem are taken to have constant values as described in the text. 
The contours correspond to masses in the range 0.1 <m< 100, with the mass increasing 
from left to right in the figure. The region in the upper left corner (above the dashed 
curve) corresponds to brown dwarfs; the region in the lower right corner corresponds to 
stars so massive that they become unstable. 

Figure 3. Empirical model for the initial mass function. Dashed curve shows the 
IMF resulting from the semi-empirical mass formula of §2 and the observed scaling laws 
which describe the distribution of effective sound speed, i.e., the distribution of initial 
conditions. The solid curve shows the Miller/Scalo fit to the observed IMF. Both curves 
are normalized to unity at m = 1 (M* = IMq). 

Figure 4. Random model for the initial mass function. Dashed curve shows the 
IMF resulting from the semi-empirical mass formula of §2 and a distribution of initial 
conditions described by a collection of n = 10 random variables. The solid curve shows 
the Miller/Scalo fit to the observed IMF. Both curves are normalized to unity at m = 1 

(M, = IMq). 

Figure 5. Random model for the IMF for different numbers n of the fundamental 
variables. The solid curve shows the Miller/Scalo fit to the observed IMF. Dashed curves 
show distributions calculated from a collection of random variables with n = 1, 2, 3, and 
5. As the value of n increases, the curves become closer to the log-normal (n 00) 
limit. All curves have been scaled so that the total width (a) of the distribution and the 
characteristic mass scale mc agree with the observed values from the Miller/Scalo IMF, 
i.e., mc = 0.095 and {a) = 1.57. In addition, all curves are normalized to unity at m = 1 
(M, = IMq). 

Figure 6. Composite distributions for the IMF using n fundamental variables with 
power-law distributions. The various curves are shown for n = 1, 3, 10 and 20. Also 
shown is the log-normal curve corresponding to the limit n — > 00. All curves have been 
scaled so that the total width (a) of the distribution and the characteristic mass scale 
mc agree with the observed values from the Miller/Scalo IMF, i.e., mc = 0.095 and 
(a) — 1.57. In addition, all curves are normalized to unity at m = 1 (M* = IMq). 

Figure 7. The effects of binary companions on the IMF. Solid curve shows the 
primary mass distribution (taken here to have the Miller/Scalo form). Dashed curve 
shows the effects of including an additional distribution of binary companions, where we 
have used the prescription of §5.3. The binary fraction JF = 0.75, the variance of the 
mass ratio distribution is a^s = 1) ^ind the mean of the mass ratio distribution is (In 0:5) 
= —1. Both curves are normalized to unity at m = 1 (M* = IMq). 
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